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Abs tra ct 


The presentation summarizes recent work in the Unsteady Aerodynamics Branch at NASA 
Langley Research Center on developing unstructured grid methods for application to steady and 
unsteady aerodynamic problems. The CAP-TSD transonic aeroelasticity code, which is based on 
the transonic small-disturbance (TSD) theory, is described first to provide background 
information to put the present work in context. The CAP-TSD code is the most fully-developed 
code for aeroelastic analysis of complete aircraft configurations at the TSD equation level and 
has been widely accepted throughout the U.S. aerospace industry. Currently, aeroelastic 
analysis capabilities are being developed at NASA Langley for the Euler and Navier-Stokes 
equations based on both structured and unstructured grids. The purpose of the presentation is to 
describe the development of unstructured grid methods which have several advantages when 
compared to methods which make use of structured grids. Unstructured grids, for example, 
easily allow the treatment of complex geometries, allow for general mesh movement for 
realistic motions and structural deformations of complete aircraft configurations which is 
important for aeroelastic analysis, and enable adaptive mesh refinement to more accurately 
resolve the physics of the flow. The presentation is therefore organized in three parts 
including: (1) steady Euler calculations for a supersonic fighter configuration to demonstrate 
the complex geometry capability; (2) unsteady Euler calculations for the supersonic fighter 
undergoing harmonic oscillations in a complete-vehicle bending mode to demonstrate the general 
mesh movement capability; and (3) vortex-dominated conical-flow calculations for highly- 
swept delta wings to demonstrate the adaptive mesh refinement capability. The basic solution 
algorithm is a multi-stage Runge-Kutta time-stepping scheme with a finite-volume spatial 
discretization based on an unstructured grid of triangles in 2D or tetrahedra in 3D. The moving 
mesh capability is a general procedure which models each edge of each triangle (2D) or 
tetrahedra (3D) with a spring. The resulting static equilibrium equations which result from a 
summation of forces are then used to move the mesh to allow it to continuously conform to the 
instantaneous position or shape of the aircraft. The adaptive mesh refinement procedure 
enriches the unstructured mesh locally to more accurately resolve the vortical flow features. 
These capabilities are described in detail along with representative results which demonstrate 
several advantages of unstructured grid methods. The presentation further discusses the 
applicability of the unstructured grid methodology to steady and unsteady aerodynamic problems 
and suggests directions for future work. 
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PRESENTATION OVERVIEW 
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CAP-TSD: COMPUTATIONAL AEROELASTICITY PROGRAM - 
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strong shock waves 




CAP-TSD INSTANTANEOUS PRESSURES ON F-16C AIRCRAFT 

DUE TO RIGID PITCHING MOTION 
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DESCRIPTION OF EULER SOLUTION ALGORITHMS 
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SURFACE TRIANGULATION 
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STEADY PRESSURE CONTOURS ON LANGLEY FIGHTER 
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OVERVIEW OF DYNAMIC MESH ALGORITHM 



Locations of points on inner boundary of grid are specified 


OVERVIEW OF DYNAMIC MESH ALGORITHM 
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ASSUMED BENDING MODE FOR LANGLEY FIGHTER 
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Mesh enrichment capability enables automatic refinement in 
regions of high flow gradients 


EFFECTS OF MESH ENRICHMENT ON CONICAL EULER 
VORTICAL FLOW SOLUTION 
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CROSS FLOW VELOCITY VECTORS FROM CONICAL 
NAVIER-STOKES SOLUTION 



a 


CONCLUDING REMARKS 



0 

> 

0 



co 

E 

0 

o 

'c 

o 

</) 

c 

CO 



CO 


0 

a> 

o 

55 

a> 


co 


CO 


o 


>» 

TJ 

k. 

0 

0 

0 

*■■■ 

0 

JZ 


0 

O) 

0 

0 

c 

o 

E 

o 

■o 

■ H 

c 

o 

0 

o> 

c 

0 


0 

k. 

X 

0 

_0 

>* 

Q. 

Q. 

■o 

3 

E 

0 

CO 

0 


o 

CO 

0 

o 

1 



■Mi 

o 

o 

0 

H— 

»*- 

0 





0 

H-* 

c 

c 


0 

o 

O 


■O£0) 
O Q. C 

£ o o 

a> o> 

E r co 

5 “ 

■U 73 Q_ 
u v. Q. 
D> 0 CO 

T3 C o 

a> 3 ■= 
3 0 ro 

O S = 

3 .2 >* 


iS >* 

3 £ 

ffl 

CO m 

o 

1- *□ 
<D 0 

3 0 
LU £ 
0 

>» C 

■O o 

0 C 

0 c 
0 

</) "O 


a - 0 
0 0 


0 


O) 

>» 

o 

■ VM 

■ MM 

c 

wmmm 

A 

o 

0 

0 

Q_ 


0 

0 

O 

Q. 

3 


0 

0 

0 

0 

E 

k. 

O 

o 


wmmm 

0 

c 

E 

0 

c 

o 

wmmm 

>» 


■o 

0 


3 

0 

O 



0 

0 

O 

c 

0 

k. 

o> 

0 

3 

T3 

0 

LLi 

4 — > 


0 

>» 

k. 

T3 

0 

0 

C 

0 

■4-^ 

o 

0 

E 

c 

0 

3 

T3 


■ 


■ 


308 


Vortex-dominated conical-flow calculations 
demonstrated adaptive mesh refinement capability 
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